function multiColumnPowersAB=computeMultiColumnABPowers(Hp, Hu, oA, cB, powersAB)
multiColumnPowersAB = zeros(Hp * oA, Hu * cB);
aux = 0;

for i=1:Hu    
    [startPointer endPointer] = getPointers(i, cB);
    
    zeroColumn = zeros(aux * oA, cB);
    gammaColumn = powersAB(1: end - aux * oA, 1:end);
    
    multiColumnPowersAB(1:end, startPointer:endPointer) = [zeroColumn; gammaColumn];
    
    aux = aux + 1;
end
end